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ABSTRACT 


A linearized, global spectral model with eight levels was used to determine 
whether the nonlinear interaction between a planetary-scale wave (wavenumber 
four) and a rapidly growing synoptic-scale wave (wavenumber seven) could 
produce a northeasterly wind, characteristic of East Asian cold surges. The 
amplitude of the synoptic-scale wave, or generic cyclone, was produced by a 
nonlinear Eady model of the atmosphere that included friction. The resulting 
nonlinear forcing was applied to either the first law of thermodynamics, the 
vorticity equation, or both. 

The thermal forcing did not produce a significant cold surge response. The 
vorticity forcing produced a respectable cold surge within 48 hours. The results of 
this study indicate the planetary-synoptic wave interaction is a possible method for 


initiating East Asian cold surges. 
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I. INTRODUCTION 


The monsoon is a three-dimensional, planetary scale wind regime that exhibits 
a strong seasonal dependence. According to the Glossary of Meteorology 
(Huschke, 1959), "the primary cause (of the monsoon) is the much greater annual 
variation of temperature over large land areas compared with neighboring ocean 
surfaces, causing an excess of pressure over the continents in the winter and a 
deficit in summer." The variation in temperature results from the position of the 
sun during each season. The shapes of the continents and their variable 
topographies produce considerable regional and temporal variability of monsoons. 

The northeasterly monsoon that occurs in East Asia during the winter is one of 
the most energetic circulations of the atmosphere. Even though the regional 
characteristics of the winter monsoon occur in East and Southeast Asia, the 
influence on other components of circulation can reach global scales (Lau and 
Chang, 1987). Because of the large scale effects on the atmosphere, the East Asian 
winter monsoon has been an area of active research. Recent areas of investigation 
include the role of the East Asian winter monsoon in midlatitude-tropical and inter- 
hemispheric interactions, monsoonal variations, and the forcing mechanisms 
responsible for small scale monsoonal variations. 

The East Asian winter monsoon 1s associated with the thermally direct Hadley 
circulation, or cell, that occurs over the area in winter. The ascending branch of 
the Hadley cell, and major monsoonal convective zone, migrates from its 
Summertime position over India to the maritime continent of Borneo/Indonesia 
(Ramage, 1971). The latent heat released in the upper levels, due to intense 


convection over the maritime continent, is transferred poleward. The cold 


Siberian high pressure system, or anticyclone, is the heat sink area for the upper- 
level poleward moving warm air. The Siberian anticyclone, іп conjunction with 
the descending branch of the Hadley cell, produces a large area of subsidence, and 
thus dominates the Southeast Asian winter (Ramage, 1971). The Himalayas block 
the southward movement of extremely cold surface air from the Siberian 
anticyclone. The only effective outflow region is to the southeast. 

Boyle and Chen (1987) documented the wintertime surface and 20 kPa wind 
fields for the period 1973 to 1984. The predominant northeasterly flow, at the 
surface, associated with the East Asian winter monsoon is shown in Fig. 1.1. The 
blacked out regions correspond to terrain heights above 1000 m. A vector length 
of 5° longitude corresponds to a ten meters per second wind, and the isotach 
contour interval is 2.5 m-s-!. The intense subtropical jet stream over Southeast 
Asia, shown in Fig. 1.2, 1s caused by the intense baroclinic zone between the warm 
tropics and frigid Siberian area. À wind vector length of 5° corresponds to a 100 
m-s-! wind, and the isotach interval is 10 m-s-!. The climatological averaged 
January 50 kPa geopotential height and wind field (Fig. 1.3) shows a dominant long 


wave trough centered over the Sea of Japan.. The contour interval is 80 m. 
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Figure 1.1. Average Wintertime Surface Wind (Boyle and Chen, 
1987) 
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Figure 1.2. Average Wintertime 20 kPa Wind (Boyle and Chen, 1987) 





Figure 1.3. Climatological January 50 kPa Height Field (Boyle and 
Chen, 1987) 


The Hadley circulation, and the effect of the Himalayas, set up the planetary 


circulation of the winter monsoon. However, Boyle and Chen (1987) indicate that 
transient synoptic-scale waves shape the final form of the Siberian High. As 
synoptic waves propagate along the longwave trough, surface cyclones and 
anticyclones develop due to the intense baroclinicity. Cyclones usually develop off 
the west coast of Japan in the area of strong upper-level positive vorticity advection 
and low-level warm air advection, and track to the northeast. Anticyclones develop 
near the southern extent of the Siberian High, due to the upper-level negative 
vorticity advection and low-level cold air advection, and track toward the 
southwest over China. When the pressure gradient between the China anticyclone 
and the cyclone off Japan tightens rapidly, significant ageostrophic motion results. 
The cross-isobaric ageostrophic flow accelerates toward lower pressure near the 
cyclone to the east and to the Intertropical Convergence Zone (ITCZ) to the south. 


The northeasterly cold wind is the "cold surge" in the winter monsoon. Cold surges 


usually reach the equatorial South China Sea in 12 to 24 hours. The enhanced 
northerly flow intensifies the tropical convection over the maritime continent, and 
thus strengthens the Hadley circulation. The cold surge ends when the midlatitude 
trough-ridge pattern moves far enough east to diminish the China anticyclone, and 
consequently the pressure gradient. Boyle and Chen (1987) emphasize that cold 
surges are caused by the interaction of synoptic and planetary waves. They 
conclude that cold surges are dynamically forced, and thus must be considered 
separate entities from the Siberian anticyclone. 

The East Asian cold surge has two stages. They are separated by a few hours to 
one day, depending on the location of the observing site (Lim and Chang, 1981; 
Chang et al., 1983). The first stage is the pressure surge. It is the leading edge of 
the air accelerating towards the equator. The pressure surge propagates with the 
speed of internal gravity waves. The second stage is a frontal passage that moves 
with advective speeds. It is defined by a sharp decrease in surface temperature and 
dew-point temperature. 

The onset of the East Asian cold surges are defined in many ways (Boyle and 
Chen, 1987). The three most common definitions are: 1) A drop in surface 
temperature at Hong Kong of five degrees Celsius, or more; 2) An increase of the 
surface pressure gradient between coastal and central China of at least 0.5 kPa; and 
3) A prevalent northerly surface flow over the South China Sea with speeds 
exceeding five meters per second(Lau and Chang, 1987). 

The cause or nature of East Asian cold surges has been examined by Lim and 
Chang (1981), Baker (1983), Bashford (1985) and Harris (1985). Lim and 
Chang (1981) used linearized shallow-water equations on an equatorial beta-plane 


to simulate the response of the tropics to a midlatitude pressure surge. They did not 


include planetary boundary layer friction or orography. Lim and Chang (1981) 
found that synoptic scale forcing in the midlatitudes produced Rossby-type waves 
that propagated into the tropics. The northeast-southwest tilt in the pressure field, 


typical of Rossby waves, and the northeasterly flow similar to cold surge events is 


shown in Fig. 1.4. 





Figure 1.4. Wind Velocity Vectors for a Barotropic Model (Lim and 
Chang, 1981) 


Baker (1983) used a global, six-layer, primitive equation model to examine the 
interaction of a midlatitude baroclinic wave with topography. Baker's results 
indicated that well developed baroclinic waves could initiate a cold surges, but the 
surges were limited, and weak. Baker concluded that other forcing mechanisms 
were required to simulate cold surges. 

Bashford (1985) used an eight-layer, spectral, primitive equation model with 
an analytical heat source to study the effects of planetary scale motion on cold 
surges. The heat source function of Bashford's baroclinic model were similar to 
that used by Lim and Chang (1981) in their barotropic model. Bashford found that 
a planetary wave (wavenumber three) with a deep thermal forcing could produce a 


cold surge response. 


Harris (1985) used the same model as Baker (1983) to study the interaction of a 
baroclinic wave (wavenumber eight) with a planetary wave (wavenumber four). 
Harris found the synoptic wave alone did not produce a cold surge, but the 
synoptic-planetary-wave interaction produced a significant cold surge response. 
Harris’ results may be tainted because the planetary wave extended unrealistically 
far south. 

This study will use a global, eight-layer, primitive equation model, and will be 
conducted in three phases. The first phase will examine the effect of a forcing with 
a single frequency derived from a forced Rossby wave, midlatitude beta-plane 
solution. The second phase will repeat Bashford's (1985) experiment with an 
impulse forcing that includes a spectrum of frequencies. This effort will attempt to 
reproduce the results of Lim and Chang (1981). The third phase will examine the 
interaction of a planetary wave (wavenumber four) with a synoptic wave 
(wavenumber seven). The synoptic wave, or generic cyclone, will be produced by 
a nonlinear Eady model of the atmosphere (Peng, 1982). The basic state or mean 
flow of the planetary wave, and the phase speed of the synoptic wave, are both equal 
to zero. Instead of combining the two waves together within the global model, as 
done by Harris (1985), this study will compute the nonlinear interactions 
analytically. The resulting wave (wavenumber three) is used in the linear, global 
model. To aid in the analysis of each experiment, the solutions from the global 


model will be projected onto the vertical normal modes. 


II. SIMPLE BETA-PLANE SOLUTION 


Since the baroclinic, primitive equation model used in this study is linear, and 


the mean or basic flow is zero, the full baroclinic atmosphere can be represented by 


the sum of vertical normal modes. The behavior of each vertical mode solution is 


similar to a barotropic shallow water system with the appropriate mean, or 


equivalent, depth (Lim and Chang, 1987). In other words, each mode of the 


baroclinic system will behave like a shallow water system. To simulate the 


variation in the Coriolis parameter, and simplify the mathematics, the system of 


shallow water equations is solved on a midlatitude beta-plane. The validity of using 


a midlatitude beta-plane is discussed by Lindzen (1967). 


The scaled shallow water equations on the midlatitude beta-plane, with a 


forcing function added to the continuity equation, are: 


О 
2 
У ф'- 10. =0 
where: 
ф' - perturbation in the geopotential height field 
V - geostrophic streamfunction 
H - mean (basic) height field or equivalent depth 
D - divergence 
V - two dimensional wind velocity vector 


C - vorticity 


fo - Coriolis parameter (constant) 

t - time 

Bo  - gradient of Coriolis parameter in the north-south direction 
y - north-south position 


- forcing function 
If the vorticity is assumed geostrophic, then ¢ = V2y, and Eq. 2.3 yields y= 
$'/fo or ó' 2 foy. When @’ is introduced into Eqs. 2.1 and 2.2, and the nonlinear 


terms are eliminated, Eqs. 2.1 and 2.2 become: 


oy _ £H 
` a 


е + Bo Ро 


When the divergence is eliminated between Eqs. 2.4 and 2.5, the 


quasigeostrophic potential vorticity equation results: 


т< 


ду 
MEN 
à gH * Bo 


Ox (2.6) 





Eq. 2.6 is solved by writing the streamfunction and forcing function as follows: 


yz W(y) e (Kx - ex) 
S-is(y)ei(Kx - a1) (2.8) 


where K is the dimensional wavenumber, and € is the frequency. The latitudinal 


structure of the forcing function is: 


| So cos? 2 y| <W 
5(у)- | 
| В УРМ (2.9) 


where So is the magnitude of the forcing, and W is half of the width of the forcing 
in the north-south direction. 
When Eqs. 2.7 and 2.8 are substituted into Eq. 2.6, the following ordinary 


differential equation results: 


d^ , gay -. 50 


dy? 9 (2.10) 
where: 
с°=- K2- fj KBo 
gH 9 (2.11) 
OT: 
-K 
ELS 
K2« o2 4 10. 
(2.12) 


The homogeneous solution to Eq. 2.10 is: 
V(y)2eztey (2 


Since the latitudinal structure is symmetrical about y = 0, Eq. 2.10 will only be 
solved for lyl € -W and y « -W. When y « -W, the negative exponent of Eq. 2.13 
is used, and when y > W the positive exponent is used. If œ? > O, the particular 


solution to Eq 2.10 for the interval lyl € W is: 


^ 
rear d =) 


e. [RE] Q.14) 
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When Eqs. 2.13 апа 2.14 are combined, the general solution (with o? » O) is: 
cos [= 

(y) - 30 d- 4 —— WI e Cei ay + ei ay) ly] < W 
2 |o? о?- [п] 

` (25155 


Wy) = Coe fay Ех (2.16) 


Since a2 > 0, a can be positive or negative. When a@ > 0, the solution to Eq. 2.16 
will have a north-south phase structure that tilts from the northeast to the southwest 
for y « -W. A northeast-southwest tilt in phase structure will propagate energy 
away from the source region. This can be shown by applying a radiation condition 
which requires that the wave have a group velocity moving away from the source. 
If & « O, the solution to Eq. 2.16 will have a phase structure that tilts from the 
northwest to the southeast, and will propagate energy toward the source region. 
For the region y > W, a < 0 will produce a northwest-southeast tilt in the phase 
structure, and propagate energy northward away from the source region. 

The two constants are determined by equating Eqs. 2.15 and 2.16, and also 
equating their derivatives, at y - -W. When Eqs. 2.15 and 2.16 are substituted into 


the real part of Eq. 2.10, the complete propagating solution is: 
192 


чу = 50 1. hw) cos (Kx -wt)- 
— | 0-7 Uu 





W 
90.1... 1 (аук (Kx + W - ot) lp «Ww 
= 
W ЕЛІ 
v -39an (aw)|-L.- — № | sn (Kx « ay- v1) y«-W 
2 о. о2-( 5: 
W (2.18) 
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If 0? < 0 and a = u7, Eq. 2.10 becomes: 


2 s(v 
E 
ау“ (2.19) 
and Eq. 2.11 becomes: 
2 
ч 
gi Oo (2.20) 
The homogeneous solution of Eq. 2.19 is: 
(у) = ену (2.21) 
and the particular solution for lyl < W 1s: 
лу 
COS E 
-..50 50 E 
ОР 2 o 
Zu a pere 
RE] (2.22) 


When y « -W, the positive exponent of Eq. 2.21 is used. 


When Eqs. 2.21 and 2.22 are combined, the general solution (with o? « 0 ог 


u2 > 0) 15: 
Ни 
COS а 
4 (у)=- 50-1 +— М | сүевуже-е) ИМ 
DANE” о 
a 
W (2.23) 
Ф (у) = Сену y«-W (2.24) 


Since [2 » 0, pt can be positive or negative. When u > 0, the solution to Eq. 2.24 
will trap energy near the source. If < 0, then the solution would grow to infinity 


for large lyl. 


The constants can be solved in the same manner as used to get Eqs. 2.17 and 
2.18. When Eqs. 2.23 and 2.24 are substituted into the real part of Eq. 2.10, the 


complete trapped solution 15: 


y — 90 гіш, 


2 Ty. 
U yu «(E 
e 1 Wcosh (uy Јн [e (Kx - wt) MES 
^ vg 
(2225) 
Y= 0sinh (pW + .___ђ___|енуст (Кх - OU y«-W 
D u? iP (ny 
w (2.26) 


The derivation above indicates that under certain conditions, or frequencies, 
the solution to the shallow water equations can either propagate energy away from 
the source region, or trap energy near the source. The experiments in Phase I will 
examine the effects of forcing the primitive equation model with a variety of 


frequencies. 
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ПІ. MODEL DESCRIPTION 


The model used in this study is a baroclinic spectral transform primitive 


equation model as described by Haltiner and Williams (1980), and used by McAtee 


(1984). The model discussion of this chapter is taken directly from McAtee. The 


model is configured to include friction, diabatic heating and a vorticity forcing. 


The specifics of how diabatic heating and vorticity are included in the model will 


be discussed in the next section. Friction is not included in this study. The basic 


equations of the model, in sigma coordinates, are as follows: 


> 








- [C f]V -K V x Бата оС 
do 

ы ү V{RTVg +6 Vio e 
о do 2 
9d. p. v.g. 0 
о до 
90 . 00 
— =- УМВ-6— + 
dt do 9 
CREE 
oo 6 
where: 
C - vorticity 
D - divergence 
T - temperature 
t - time 
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+kVxXF+A 





(° 


(258) 


(5 


(3.4) 


(39559 


- potential temperature 

- Surface pressure 
horizontal velocity vector 
- geopotential height 


- gas constant 


с И секо ие ет ар 


М - specific heat at constant pressure 





f - Coriolis parameter 
G - vertical coordinate (O = E EE 
II - PT 
| | ас 
б - vertical velocity (6 — Eme 
а - nA 
- pressure 
Pr - pressure at the top of the model 
Е - frictional force 
- diabatic heat forcing function 
A - vorticity forcing function 


The continuity equation (3.3) is rewritten by integrating with respect to sigma. 
and applying the boundary conditions 6 (0) 2 6(1) = 0.0. Thus the integral of 
Eq. 3.3 may be written: 


(520) 


15 


The vertical velocity, 6, 1s derived diagnostically by substituting Eq. 3.6 into 


Eq. 3.3, integrating in the vertical, and using 6 (O) = 0. 
5-5 | (D+ Glo 
0 (3.7) 


The first law of thermodynamics, Eq. 3.4, is written: 


9T. v. VT-o* одат) кт +Q 
О 


99 ууа 
o Ot 








(3.8) 


where x equals R/Cp, To apply semi-implicit differencing, it is necessary to divide 


the temperature into parts as follows: 


T=T*(o)+T o). @ t) (3.9) 
where: 

I: - appropriately averaged temperature 

T - perturbation temperature 

À - longitude 

© - latitude 


The basic equations are conveniently written in spherical coordinates by 
defining the following operator: 


26 ل 


—— њи њи — M — — 


с соѕ2ф Әх Ф Әф (3.10) 


Using Eqs. 3.9 and 3.10, the basic equations can be written as follows: 


АВ 

ot (3.11) 
DUE SA TOR. +RT*q} 

à (3.12) 
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-a (UT, ID o5 9——— 


XTc") | e r(o. G-D)«Q 
do 


(3.13) 


(3.14) 


(3.15) 


where: 


RT COS 
A=(C+ нда соёр d. rudi 





pU V 
2со$-Ф 


_ №005 Ф 








и - zonal component of the velocity vector (V) 

у - meridional component of the velocity vector (V) 

Т - mean radius of the Earth 

Eqs. 3.10 through 3.15 are the basic equations used in the model. These equations 


are represented spectrally in the horizontal and finite differenced in the vertical. 
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A. VERTICAL STRUCTURE 


The vertical structure of the model follows the development given by Arakawa 


and Suarez (1983). The variables are staggered in o so that G, D, U, V and T are 


carried at the mid-point of each layer, where 6 =6,. The variable 6 is carried at the 


top and bottom of each layer, where o= ok. The vertical structure is illustrated in 


Fig. 3.1. The finite difference form of Eqs. 3.10 through 3.15 are: 


02-a (AB) + A 





т OB AT. Ук +RT*q +E] 
A -(G+D} 
Tk =o, (UT.VT]T +T,(G,-G-D)- 
Len т OAK (Th ысты |+0 








Qk - Qk+1 = ОЛБО - н: 4 ВкТкет 
Pk Рк+1 





* K 
QLM = Qs + CPTIM | 
PLM 


Oj.1 2 Op (G4 D)- У (с j + Dil^o; 








ЈЕ! 
where: 
P. = (oxa) 
Полкот“) 
Р= 1+к 
Ok+1 - Ok 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


(3.2 


(3 





| ЕТ'усо$2ф | 1 94 (6% 
A= (++ Bhs br ViVi) [RT Ce ст 
HOO] 


| cos Q 
Oy Сы E Ux) + oi (Ux а Чл) Е сте] | 8 јакој 
6 г (оф 00 


V, VII 
Ш Пп 


С 


D=V-V; 


Бер 


E - 
k Pk+1 - Pk 


_ Рр - Ра = | L A 
"E^ Pia - Pi 
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Figure 3.1. Vertical Structure of the Primitive Equation Model 

Eqs. 3.16 through 3.22 are written in matrix form, so the terms on the right 
hand sides contain all the terms that are to be evaluated explicitly, and the terms on 
the left hand sides contain those terms that are to be evaluated implicitly. Eq. 3.20 


is combined with an integrated finite difference form of Eq. 3.15 to obtain: 


ф=|С|Т+ф; (3.23) 


where ICI is a square matrix and the other quantities are column vectors. The finite 


difference form of the surface pressure tendency equation, Eq. 3.14, is : 


a- > (Ок + Dy Ack 
К (3.24) 


which can be written 1n matrix form as: 
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9d... NTIG.D) | 
д (3.25) 


where NT is the transpose of a constant column vector. Similarly, Eq. 3.22 is 
written: 


6=ZAG+D) (3.26) 


The next to last term in Eq. 3.19 is: 





а. ав Fk Кал - тА Pena] 


jr k-1 O 


For the purpose of semi-implicit formulation, the temperature is separated 


according to Eq. 3.9. The mean part of that term is written: 











Sabi Pia - mea [i- Beni ]=IM](G+D) 

I Ру (9:20) 

Ед$. 3.17, 3.19, 3.23 and 3.25 are now written: 

ШЕ) кг -а(в,А)-У-Е-Кр 

д (3.29) 

4 + ND =-NIG 

д (3.30) 
Т 

IT Мор = Kr 

à (6300) 
>’ =|C|T (3.32) 


where IQ] = IMI + KT*NT and 6' = ф —- ós. The variables Kp and Kr represent terms 
that are not explicitly separated out. 
The semi-implicit time differencing is achieved by evaluating the terms on the 


left hand sides of Eqs. 3.29 through 3.31 implicitly. The remaining terms and 
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Eq. 3.32 are evaluated explicitly using leapfrog differencing. The difference 


equations are written: 


2 
Dij F AtV (CT F | R Tee 


ро ДК C| T, 4 | RI T*q, i) - 2At(Kp), 


(3.33) 
Т + АО рь = Tn - А0010 - 2At (Kr) (3.34) 
Ф) + АКТ = да] AND ANGE (3.35) 


Now, the following equation for D is found by substituting Eqs. 3.34 and 3.35 into 
Ра 95 


В. ы  B,Ds1 * 2A (Kp AVET | R| T*q, 4 | C] (Kz),)- 
IR| T*NTG, (3.36) 


where the matrix operator B is: 


B. -AX CIIQIE RITNT)V" +|1 (3.37) 


and IIl is the identity matrix. 


B. SPECTRAL FORMULATION 
The equations 3.10 through 3.15 are represented spectrally in the horizontal. 


The variables are represented as follows: 


E EE 
C(A.o.o:]- Y, Y, Mosom- Y) У сту 
т--) пет) та-) пет) (3.38) 


where C is some variable, and: 


СШЦ)-(-1) СЕЕ 


where m is the zonal wavenumber, n is the meridional index, and n —|ml gives the 
number of zeros between the poles (-1 < sin @ < 1) of the associated Legendre 
function. Triangular truncation is used in this study, with the truncation limit, J, 
equal to three. The non dimensional zonal coordinate index, A, equals (5 – 1)/2 
where 1 < s < 16. Note that the separation is such that the coefficients C5 are 
functions of time and the vertical and spherical harmonic YT are horizontal 
functions of space. The normalization and orthogonal properties of YF allow the 


coefficients to be obtained as follows: 


жо | 

CF =(C, Y= | | СупафдА 

0041 (3.39) 
The nonlinear terms are computed using the transform method following Haltiner 
and Williams (1980). The longitudinal direction is treated with a Fast Fourier 
Transform and the latitudinal direction uses Gaussian Quadrature. The number of 
latitudes, N, and longitudes, M, satisfy: N> 3J/2+ 1 and M > 3(J -1) +1. For this 
study, N = 60 and M = 48. 


IV. FORCING FUNCTIONS 


Chapter III of this paper describes the global spectral model. The vorticity 
equation (3.1) and the thermodynamic equation (3.4) use forcing functions, A and 


Q, respectively. This chapter describes the forcing functions. 


A. PHASE I - ROSSBY WAVE FORCING 
The thermal and vorticity forcing of Phase I follows the development of 


Chapter II. The forcing function is: 
F=(Fa)(Fs)(Fo)(Fr) (4.1) 


where Fa is the amplitude of either the thermal or vorticity forcing. The vertical 


Structure (Fs), shown in Fig. 4.1, 1s: 


_ 2 2 
| (1.0 Len OLO 026m. 
sinh“ 1.0 
Бе = 
tanh У Он (ОН G S Om 
tanh 1.0 (4.2) 


where Fao 1s the amplitude of either forcing at the surface (for this study Fag = 0), 


and: 
O07 O - Omax 

1.0 - Стах (4.3) 
С = О - Omin 

Стах - Omin (4.4) 
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Figure 4.1. Vertical Structure of Forced Rossby Wave Solution 


The meridional structure (Fo) is: 


Era ex Е Үз) 
YN- Ys (4.5) 

where @ is the latitude, and YN — Ys is the north-south period of forcing. 
The development in Chapter II uses a cos? latitudinal function, where the actual 
latitudinal forcing of the model is sin. The difference is insignificant. The forcing 


is equal to zero in the Southern Hemisphere, and outside of the latitudinal limits in 


the Northern Hemisphere (Fig. 4.2). 
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Figure 4.2. Latitudinal Structure of Forced Rossby Wave Solution 
The amplitude, vertical and latitudinal components of the forcing function 
discussed so far are similar to Bashford (1985). Both the temporal and longitudinal 
structure of the forcing functions are included in the function Fg shown below: 


эы a 


А ХЕ Хи 


(4.6) 


where Xw — Xg is the East - West period of forcing, and A is the longitude. The 


frequency (0) is given by Eq. 2.11. 


B. PHASE If - IMPULSE FORCING 

Bashford (1985) used an impulse forcing in the thermal equation similar to 
Lim and Chang (1981). In this study the impulse is also added to the vorticity 
equation. Since this impulse forcing includes a spectrum frequencies, it is more 
realistic than the single-frequency forcing described in Phase I. The impulse 


forcing is given by: 
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F- (FAJ(Fs)(Fo)(FA)(Fr) (4.7) 


where Fa, Fs, and Fg are the same as Phase I. The longitudinal part (F5) is the 
same as Phase I, but it is separate from the temporal function. The impulse is 


generated by the temporal function (Fr): 


Fr=- et 
213 (4.8) 


where t is time and Tis the time scale. The impulse peaks at t = 2T and then decays 
exponentially to zero. The solutions in this study use Tt = .5 day, so the peak of the 


impulse occurs at one day (Fig.4.3). 
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Figure 4.3. Temporal Structure of the Impulse Forcing 
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C. PHASE III - GENERIC CYCLONE FORCING 
1. Thermal Forcing 
The thermal forcing from the generic cyclone represents the nonlinear 
advection of the synoptic temperature field by planetary waves, and vice versa. 


The nonlinear advection terms are given below: 


Q=- Vr -VTs - VWs: VTL (4.9) 
where: 

VL - two dimensional planetary wind vector 

Vs - two dimensional synoptic wind vector 

T. - planetary scale temperature field 

Ts - Synoptic scale temperature field 


Eq. 4.9 is written in scalar form as follows: 
OTs OTs ӨЛІП! OTL 

Qz-U S. Vi S. Ug L. v E 
Ох ду Ох ду (4,10) 
Тһе quasistationary geopotential height field for the planetary wave is: 

OL = A(P,t)W(y)cos (Kx) (4.11) 


where A(P,t) 1s the amplitude factor for the planetary wave (constant for this 
Study), and K_ 1s the planetary wave number (for this study KL = 4). The latitudinal 


Structure of the planetary wave is: 





Wty) =n {42 Yss] 
Үрм - Yps (4.12) 
where Yps — Ypnis the north-south period of forcing for the planetary wave. 
The geopotential height field for the synoptic wave is: 
фа = Муја Ров (Кох - vt) - b(P,t)sn (Ksx - vt) (4.13) 


where a(P,t) and b(P,t) are the vertical amplitudes of the synoptic scale wave. The 
functions, a(P,t) and b(P,t), represent a baroclinic wave generated from an Eady 
model, to be discussed later. The value Kg is the synoptic wavenumber (for this 


study Ks = 7), and v is the frequency. The latitudinal structure of the synoptic-scale 





wave 15: 
i - Y 
N(y)=sin тіу s 
NGN = Yss (4.14) 
where Ysn — Yssis the north-south period of forcing for the synoptic wave. 
The ideal gas law: 
рК (4.15) 
and the hydrostatic equation: 
ОР 
=e. PE 
97 (4.16) 
are combined to yield: 
д 
T< =- POS _ Ра (Кох - мо + Фа (Кох - vt} 
ROP К (ОР ОР (4.17) 
апа: 
О 
T=- L1 POA Wes (Ку х) 
К oP К ӘР (4.18) 


The planetary wave amplitude function, A(P,t), is constant for this study 


so TL = 0 and Eq. 4.10 simplifies to: 


Т 10 
О= 11915 у 915 


Ох ду (4.19) 


The planetary and synoptic vorticity fields are: 


DO 


и 
(4.20) 


9 М 


5-19; = AE. кім [acos (Ksx - vt) + bsin (Ksx - vt) 





(4.21) 
The planetary and synoptic zonal wind fields are: 
ТЕКШЕ AW cos (Kr x) 
foy f а (4.22) 
0 “7 
Us =- 1995. и “Pa. (Кох - vt) +bsin (Ksx - vt) 
{ду 2 (4.23) 
The planetary and synoptic meridional wind fields are: 
v= Lok = AWK; sin (Kix) 
fox f (4.24) 
_ 1995 _] Ub 
Vs = + = NKs|[-asin(Ksx - vt} +beos (Ksx - №1] 
f ox f (4.25) 


When Eqs. 4.17, 4.22 and 4.24 are substituted into Eq. 4.19, and 
simplified, the result is: 
Q=-F A РА (макоо кух) (Ksx - vt) + 


pr (Ki x Jeos (Kex - vt)+ 


“кэп ОООО О m 


gr o | (Кух јат (Кох - v J 


dy ФР (4.26) 
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Using the trigonometric identities: 


sin (a - B ) - sino cos B - cos a. sin 
cos(a -B ) 2 cos o. cos D sin o: sin B 


Eq. 4.26 becomes: 


о--РА ОМ, „оа 
Rf L oy 


ОУ „о 


57 Ме Кссов (Ку х (00 (Кох cos (vt) + sin (Кох ап (уп) + 


ду др 


s KS хўох (Ksx (os (vt) + sin (Ksx)sin (vt) + 





WA ki Ban (Krasin (Ksxjeos (vi)-cos(Ksxhin (vi) | 


Using the trigonometric identities: 


sin o.00s B= U{ sin (cr +B) + sin (cr-B ) 
cos a cos B = (cos (o.+8) +cos(o.-8 ) 
sn asin В = (cos (o-8 )-cos (0+8) 


Eq. 4.29 becomes: 


B 


= Na Kscos (Ki x(sn (Kgx)eos (vt) - cos (Ksx)sn (vi)) 


(4.27) 


(4.28) 


(4.29) 


(4.30) 


(4.31) 


(4.32) 


Ое | SNK [суй (Ks + Ki) x) + sin (Ks - Kj )x))- 


ок eco (üt ДЕ 


мак, E (vikcos (Ks -- Ki x) coss - Ki)x) 


sin (vifsin ((Ks + Ky )x) +sin (Ks - Ku) | + 
we ки aE (vt\sin (Ky, + Ks)x)-sin (Ks - Кук) + 
sin (vifcos ((Ks - Kx 5x) -cos (Ks - кА) | + 


к EIS Ki x) -eos (Ks Ki )x)) 


sin (vtXsin (Ky + Ks)x)-sin ((Ks - KJ) | | 
(4.33) 


Eq. 4.33 describes the nonlinear interaction between the planetary-scale wave 
(wavenumber four) and the synoptic-scale wave (wavenumber seven). The 
resulting waves are wavenumber three (Ks — KL) and wavenumber 11 (Ks + Kr). 
Since synoptic-scale waves do not propagate into the tropics (Harris, 1985), 
wavenumber 11 is ignored in Eq. 4.33. The resulting equation, (4.34), is the 


thermal forcing function, Q: 


о--1БА | T NAKS E (visin (Ks - Ki )x)- sin (vtjcos (Ks - к + 


к Е (vicos (Ks - Kx x) sin (vilsin (Ks - ко | + 


к [ cos (vt)sin (Ks - Kx) sin (vijeos (Ks - ко | + 


yg, COS (Усов ((Ks - К; )х) t sin (уып (Ks - Ku) | | 


oP 
oy (4.34) 


2. Vorticity Forcing 
The vorticity forcing function represents the nonlinear advection of 
vorticity by both the planetary-scale and synoptic-scale waves. The nonlinear 


advection terms are: 


A =- Vi -VÚs - Ve Vig (4.35) 


Eq. 4.35 is written in scalar form as follows: 


5077025: 022,90 
De. Vy - Us - Vg 
Ox ду Ох ду (4.36) 


When Eqs. 4.20, 4.21, 4.22, 4.23, 4.24, and 4.25 are substituted into 
Eq. 4.36, the result is: 


oe 


aN -K ык asin (Ksx - vt) + bcos (Kgx -у1))+ 


م 
ду dy?‏ { 2{ 
WK aN KEN |g‏ 
әу 79‏ 


ƏN.. Jlo W 
acl? ask 





sn (Kj xYacos (Ksx - vt) + bsin (Ksx -vt))- 





sin (Kj xYacos (Ksx - vt) + bsin (Ksx -vt))- 





S aw : к20У (Ki xY- asin (Ksx - vt) + Бсо5 (Кох ZI 
33 ey 


NK 





(4.37) 


When like terms are grouped together and the trigonometric identities in 


Eqs. 4.27 and 4.28 are used to simplify Eq. 4.37, A becomes: 
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E GP 











OW QW 20W 
val [RE ú хау ка | 


Б (Kix)(- asin (Kgx cos (wt) + acos (Ksx)sin (vt) + 


bcos {K sx kos (vt) + bsin (K5x)}sin (ғ) | ES 


ON 629“ ӘМ IW 
ШЕ | а к 


Б (Ку х)(асоѕ (Ksx cos (vt) - asin (Ksx sin (vt) 





iss) ts nt) | | 


(4.38) 


Using the trigonometric identities of Eqs. 4.30, 4.31 and 4.32, Eq. 4.38 
becomes: 
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D 





^) 
- KEN |- 


дуг 





- NK aw ge on 
E ду 








Ы. 


[ asin ((Ks + К1/х)-авл(К- Кі/х/сов(у1)-- 


| 


(acos (Ks + K1 )x) + acos ({(Ks - K1)x)}sin (vt) + 
(teos ((Ks + K1 )x) + bcos ((Ks - Ki )x))eos (vt) ^ 


(bsin ((Ks + K1)x) + bsin ((Ks - K1)x)}sin | + 


Бы aN QN | ON OW a ] 
ày a] Oy lay 


| n (Ко + Кијх) - asin (Ks - Ku)x)cos (vt) + 








(acos ((Ks - Ку Јх) - асо5 (Ко ч- Kr)x)JJsin (vt) 


(bcos (Ks - Кі/х)-Әсов (Кө > Кі/х/ров(м() 


Eq. 4.39 describes the nonlinear interaction between the planetary wave 


(bsin (Ks -- Ki )x)- bsin (Ks - Ki )x) 
(4.39) 


(wavenumber four), and the synoptic wave (wavenumber seven). The resulting 
shortwave, wavenumber 11, is ignored. The final form of the vorticity forcing 


function is: 
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aN N 
~- KgN 








-NK АД ү 
дуз ду 





[> | 


| asin (Ks - Kr )xjcos (vt) * acos (Ks - Kr)x)sin (vt) + 


bcos ((Ks - Ki.)x)cos (vt) - bsin ((Ks - Ki xen (vt) | + 


(4% | 


| asin ((Ks - Kr xjeos(vt) - acos ((Ks - Ki x)sn (vt) + 





к ju? e KiW 





bcos ((Ks - K)x}cos (vt) + bsin ((Ks - Ki )x)sin J | 
(4.40) 
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V. DESCRIPTION OF GENERIC CYCLONE MODEL 


The synoptic-scale geopotential height field is given by Eq. 4.13, where the 
vertical amplitudes, a(P,t) and b(P,t), are obtained from the growth of a generic 
cyclone. The intensity of the generic cyclone follows the development of 
midlatitude baroclinic disturbances in the atmosphere, using an Eady (1949) model. 
Eady assumed the atmosphere was a frictionless, continuously stratified fluid 
whose motion is adiabatic, hydrostatic and quasi-geostrophic. Eady’s original 
study neglected the nonlinear wave interactions. 

Peng (1982) applied spectral methods to the Eady (1949) model, but included 
frictional dissipation and nonlinear wave interactions. Peng found the behavior of 
the baroclinic wave depended on the stratification, S, frictional dissipation, Y, 
supercriucality of the vertical wind shear, A, and the fundamental zonal wave 
number, k. The Peng (1982) model of a baroclinic wave is used in this study to 
generate the generic cyclone. 

The atmosphere is considered to be an infinite channeled Boussinesq fluid with 
a constant Brunt- Vaisala frequency, Ng. The top and bottom boundaries are rigid 
walls, separated by a distance, D. The lateral boundaries are also rigid walls, 
separated by a distance, L. The basic state is a zonal flow, U, with a constant 
vertical wind shear, à. The quasi-geostrophic potential vorticity equation that 


describes the motion 1s: 


aes 
CN EON C T 
nu га (wq) 





(5.1) 


where: 
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_дхду дудх (5.2) 
and the potential vorticity (q) is: 


2 2 2 
-9у 9у 19У 


дх2 ду: 5 922 (5.3) 


with y being the disturbance streamfunction. The vertical variable (z) ranges from 
minus 1/2 at the bottom of the model, to 1/2 at the top. 


The basic state stratification parameter, S, of the fluid is constant,and given by: 


;- ND? 
612 (5.4) 
where: 
Ng - Brunt- Vaisala frequency 
D - vertical depth of the fluid 
fo - Coriolis parameter 
j: - width of the channel 


The frictional dissipation, Yy, is given by: 
Ro (055) 


where E, is the Ekman number, and Rg is the Rossby number. 
The supercriticality of flow, A, measures the vertical wind shear in excess of 


the critical vertical wind shear required for growth of the baroclinic wave, and is: 


А= А-А, (5.6) 


where: 
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_ 4y u?tanhż u 
k [(}1 - tanh )(1 - Hanh i J (5.7) 


represents the marginal shear required for instability by linear theory, and: 
1 
ЧА? 
=k +2} 


where k is the fundamental zonal wave number. The marginal stability curve, as a 


(5.8) 


function of [L for S = .0628 is given in Fig. 5.1. 
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Figure 5.1. Marginal Stability Curve of Eady Model, Including 
Friction 


The disturbance potential vorticity (q) is equal to zero throughout the lifetime 


of the cyclone in this model. Eq. 5.1 simplifies to: 
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2 2 2 
U OU Јом 











Әх? 2 S 97? (5.9) 
The nondimensional boundary conditions for the model are: 
Wo at y =0,1 
Ох (5.10) 
2. 

Пра у =, 

дуд (5.11) 
зе ге! (5.12) 


where V is the zonal average streamfunction. 
Eqs 5.9, 5.10, 5.11 and 5.12 form a closed system of equations. A Fourier 
series, Eq. 5.13, whose y structure satisfies Eqs. 5.10 and 5.11, is used to represent 


the horizontal streamfunction field: 


у= У, У, (Солей + Сет) хіп плу + У, Съсоѕ пту | 
m n n (5.13) 


where m and n are the number of zonal and meridional modes, respectively. The 
complex amplitude coefficient (Cinn) of mode (m,n), and the complex conjugate 
(C* n) Of Cmn, are both functions of time. Since only one baroclinic wave is used 
in this study, Eq. 5.13 simplifies to: 

y=(C, e + С\үе-®х) sin Ty + Cocos Ty (5.14) 


The amplitude coefficients, Cj;, C*,;, апа Со are determined by substituting 
Eq 5.14 into Eqs. 5.9, 5.10, 5.11 and 5.12. The system of equations is integrated in 


time using a fourth order Runge-Kutta method (Peng, 1982). The vertical 


4] 


amplitude coefficients, a(P,t) and b(P,t), in Eq. 4.13 are related to the amplitude 


coefficients of Peng by: 
Cı1 = acosh (uz) + b sinh (pz) (5.15) 


The variables a and b are determined by solving Eq. 5.15 atz 2 € 1/2. 

The nondimensional coefficients a and b are made dimensional by multiplying 
the them by 4 000 000 m (the characteristic length scale of the Eady model), 40 
m-s-1 (the characteristic velocity scale), and .000 1 s-! (the Coriolis parameter at 


45° №). 


VI. ANALYSIS OF RESULTS 


A. NORMAL MODE ANALYSIS 

The solutions of all three phases of this study will be analyzed using the vertical 
normal modes of the primitive equation model. Gill (1982) provides a good 
general discussion of normal mode analysis, and Lim and Chang (1987) derive the 
vertical modes of a shallow water equation model solved on an equatorial beta- 
plane. 

The vertical normal modes of the primitive equation model used in this study 
are displayed in Figs. 6.1 to 6.8. There are eight vertical modes, consistent with the 
eight levels of the model. The first mode does not have a zero-crossing. and is 
considered a barotropic mode. The second mode has one zero-crossing. Since the 
amplitude of the second mode changes sign in the vertical, it is considered the 
baroclinic mode (Lim and Chang, 1987). The higher modes have successively 
more zero-crossings. The solutions to the primitive equation model at each level 


are projected onto these vertical profiles and summed to produce the modal output. 


The equivalent depth (H) of each vertical mode is provided in Table 6.1. 
TABLE 6.1. EQUIVALENT DEPTH OF VERTICAL MODES 


Vertical Mode Equivalent Depth, m 
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Figure 6.1. Vertical Profile of Normal Mode 1 
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Figure 6.2. Vertical Profile of Normal Mode 2 
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Figure 6.3. Vertical Profile of Normal Mode 3 
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Figure 6.4. Vertical Profile of Normal Mode 4 
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Figure 6.5. Vertical Profile of Normal Mode 5 
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Figure 6.6. Vertical Profile of Normal Mode 6 
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Figure 6.7. Vertical Profile of Normal Mode 7 
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Figure 6.8. Vertical Profile of Normal Mode 8 
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B. PHASE I RESULTS 

The forcing functions of Phase I will use the development of forced Rossby 
waves from Chapter II. In each experiment the forcing is applied between 24° N 
апа 36° М, and the model is integrated forward in time for 240 hours. The 
frequency (€) and corresponding period of the temporal component of the forcing 
function (Eq. 4.6) are given in Table 6.2. Recall from Chapter IJ, that the sign of 
the variable a? determines whether the solution to each vertical mode will 
propagate away from the source region. The variable o2, as a function of 
equivalent depth and frequency, is given in Tables 6.3a and 6.3b. The o? values are 
determined by first setting &2K 23.797 x 10-7 m-l. This value of K 
corresponds to wavenumber three at 30° N. The frequency is calculated for each 
equivalent depth by using this value of K in Eq. 2.12. The variable a? is then 
calculated for each frequency and equivalent depth using Eq. 2.11. Both fo and Bo 
of Eqs. 2.11 and 2.12 are calculated at 30° N. When a2 > 0, the solution for that 
specific equivalent depth and frequency will propagate energy away from the 
source region. According to theory, only the barotropic mode (H = 7083.1 m) will 
propagate energy away from the source when the frequency (@) equals 
-206.30 x 10-7 s-1. As the frequency (Qo) decreases, the higher modes have 


propagating solutions. 
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TABLE 6.2. ROSSBY WAVE FREQUENCY AND PERIOD OF 
VERTICAL MODES 


Frequency (a 10- 


Period, days 





TABLE 6.3A. ALPHA SQUARED (M-! X 10-12) FOR HIGH 
FREQUENCY MODES 


Frequency (@), s-! x 10-7 
-206.30  -26.13  -8.04 -- 3.88 


7083.1 2.60 9.14 19.18 
209.1 0.14 6.62 16.66 
59.7 634 0.14 10.18 

Depth 28.4 -16.38 -9.90 0.14 
(H), 16.1 -30.86 -2438  -14.34 
m 9.8 —5236 —588 4 

6.1 -85.29 -78.81 68.77 

3.8 -140.00-133.50 -123.50 


TABLE 6.3B. ALPHA SQUARED (M-! X 10-12) FOR LOW 
FREQUENCY MODES 


= 20090 


Equivalent 





Frequency (0), s-! x 10 -7 


-1.36 -0.85 -0.53 *206.30 


Equivalent 


Depth 
(Н), 


* changed the sign of @ associated with mode 1 


7083.1 
209.1 
39.7 
28.4 
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88.10 
85.58 
79.10 
69.06 
54.57 
33.08 
0.14 
-54,59 


142.80 
140.30 
133.80 
123.80 
109.30 
87.81 
54,87 


-0.59 
—3.10 
-9.58 
-19,62 
—34.11 
—55.60 
-88.54 





0.14 -143.30 


Phase I consists of seven experiments. Each experiment will help determine 
how the primitive equation model reacts to the analytical forcings derived in 
Chapter II. The seven experiments are summarized in Table 6.4. The first three 
experiments will determine the effects of a thermal forcing with different 
equivalent heights (H) or frequencies (o). Experiments 4, 5 and 6 will determine 
the effects of a vorticity forcing with different frequencies. Since the sign of @ in 


Experiment 5 is reversed, 02 is negative, and all Rossby waves should be trapped 


near the source. The variable & in Experiment 7 is a singular point in Eq. 2.15. 


TABLE 6.4. PHASE I EXPERIMENTS 


Experiment Thermal Forcing Vorticity Forcing Equivalent Depth 





* sign of o changed 
жо = лу М 


1. Experiment 1 
In this experiment the forcing function is only added to the thermal 
equation (Eq. 3.4). The mean or equivalent height (H) in Eq. 2.12 is 7083.1 m. 
which is exactly equal to the scale height of the first vertical mode (the barotropic 
mode). The wave frequency (o) is 2206.30 x 10-7 s-! and the period is 3.5 days. 
The wind velocity vectors for the lowest level (o = .938) are displayed in Figs. 6.9 


to 6.18. The scale of the wind vectors is shown in meters per second at the bottom 


of each figure. 
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Figure 6.10. Wind Vectors at 48 Hours for Phase I, Experiment 1 
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Wind Vectors at 72 Hours for Phase I, Experiment 1 


Figure 6.11. 
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Figure 6.12. Wind Vectors at 96 Hours for Phase I, Experiment 1 
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6.13. Wind Vectors at 120 Hours for Phase I, Exper 
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Figure 6.14. Wind Vectors at 144 Hours for Phase I, Experiment 1 
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Wind Vectors at 168 Hours for Phase I, Experiment 1 


Figure 6.15. 
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Wind Vectors at 192 Hours for Phase I, Experiment 1 


Figure 6.16. 
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Wind Vectors at 216 Hours for Phase I, Experiment 1 


Figure 6.17. 





Wind Vectors at 240 Hours for Phase I, Experiment 1 


6.18. 


Figure 
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In the first 24 hours it is evident that wavenumber three develops, since 
there is one wave in each 120° sector. The circulation is confined to the latitudes of 
the forcing. By 48 hours the overall magnitude of the vectors has increased, but 
there still is not significant propagation into the tropics. At the 72 hour point the 
magnitude of the wind vectors has decreased, and the circulation has pushed 
southward toward the equator. The wave patterns move toward the west, but the 
expected northeast-southwest tilt of Rossby waves is not present. By 96 hours the 
overall magnitude of the wind vectors is at the minimum, and there 1s weak 
circulation extending into the Southern Hemisphere. After 120 hours the overall 
magnitude of the wind vectors has increased to the value at 24 hours, and the weak 
circulation in the tropics is not evident. 

A similar pattern of changes occur at 144, 168, 192, 216 and 240 hours. 
The three-and-a-half day oscillation, or beat, of the overall velocity vector 
magnitude is consistent with the three-and-a-half day period corresponding to the 
frequency of the waves in this experiment; however, the beat frequency is not 
consistent with the beta-plane theory. The oscillation in magnitude indicates the 
wave pattern has not stabilized, and the time variation is not purely sinusoidal, as 
assumed in the theory. 

The circulation patterns in Figs. 6.9 to 6.18 do not propagate very far 
north. The trapping of waves to the north may be due to an increase in the Coriolis 
parameter on the beta-plane. As the Coriolis parameter (f) increases in Eq. 2.11, 
a? becomes more negative; consequently, the waves are trapped near the source 
latitudes. 

The wind vectors for the first vertical mode at 24 and 240 hours are 


displayed in Figs. 6.19 and 6.20. At the 24-hour point, there is significant 
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circulation in the tropics, but the magnitude is only one meter per second (one fifth 
of the complete field shown in Fig. 6.9). This observation is consistent with theory, 
in that it is difficult for a heating or mass source to force the barotropic mode (Lim 
and Chang, 1987). After 240 hours, the overall magnitude of the wind vectors has 
doubled. There is significant circulation south of the equator, but the desired 


northeast-southwest tilt of propagating Rossby waves is not evident. 
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Figure 6.19. Wind Vectors at 24 Hours for Phase I, Experiment 1, 
Mode 1 
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Figure 6.20. Wind Vectors at 240 Hours for Phase I, Experiment 1, 
Mode 1 


The pattern of the first baroclinic mode at 24 hours (Fig. 6.21) is similar to 
the barotropic mode at 24 hours, but the overall magnitude of the wind vectors is 
larger. This result is similar to Lim and Chang (1987). Lim and Chang found that 
a heating source is more effective at forcing the higher order modes. The 
circulation pattern and overall magnitude of the wind vectors of the first baroclinic 
mode, after 240 hours (Fig. 6.22), 15 similar to the barotropic mode at 240 hours, 
except the modes are 180° out of phase. The overall magnitude of the wind vectors 
of the higher modes (not shown) are less than or equal to the magnitude of the 
barotropic mode at 24 hours. According to theory, the barotropic mode is the only 
mode with a propagating solution (when H = 7083.1 m). The results of this 


experiment indicate all modes can propagate to some extent. 
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Wind Vectors at 24 Hours for Phase I, Experiment 1, 


Figure 6.21. 
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Wind Vectors at 240 Hours for Phase I, Experiment 1, 


Figure 6.22. 


Mode 2 
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2. Experiment 2 

In this experiment the forcing function is again added to the thermal 
equation, (Eq. 3.4) but the equivalent depth (H) is 209.09 m. This value of H 
corresponds to the first baroclinic mode, which 1s mode two. The wave frequency 
(w) is 26.13 x 10-’ s-! and the period is 27.8 days. The wind vector fields for the 
lowest level (o = .938) at 24 and 240 hours are displayed in Figs. 6.23 and 6.24. 
The patterns are not significantly different than Experiment 1. The overall 
amplitude of the wind vectors is much larger than for Experiment 1. The beat 
oscillation of the overall vector magnitude is not present, since the period is 27.8 
days. The westward propagation 1s slower than in Experiment 1, also due to the 
longer period. The overall magnitude of the wind vectors of the barotropic mode 
at 24 hours (not shown) 1s one meter per second, whereas the overall magnitude of 
the baroclinic wind vectors (not shown) 15 ten meters per second. The magnitudes 
of the wind vectors of the higher modes are also on the order of the barotropic 
mode at 24 hours. All the modes at later times show similar patterns and 


magnitudes. 
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Figure 6.23. Wind Vectors at 24 Hours for Phase I, Experimen 
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Figure 
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3. Experiment 3 

The forcing function of this experiment is added to the thermal equation 
(Eq. 3.4), and uses an equivalent depth (H) of 59.748 m. The frequency (0) is 
-8.04 x 10-7 s-! and the period is 90.4 days. The results are similar to Experiments 
] and 2, except the overall magnitude of the velocity vectors is greater and the 
westward propagation of the circulation patterns is slower. The magnitude of the 
wind vectors of the barotropic mode is again much less than the baroclinic mode. 

4. Experiment 4 

In this experiment the forcing function is only added to the vorticity 
equation (Eq. 3.1). The equivalent height (H) in Eq. 2.12 is 7083.1 m, which is 
the height of the barotropic mode. The frequency (œ) is 206.30 x 10-7 s-! and the 
period is 3.5 days. The wind vector fields for the lowest level (o = .938) are 
displayed in Figs. 6.25 to 6.34. The circulation patterns are substantially different 
than when the source is added to the thermal equations in Experiments 1 through 3. 
After 24 hours there is significant circulation moving toward the tropics. By 48 
hours the flow has crossed the equator and has the northeast-southwest tilt 
consistent with propagating Rossby waves. The 72 through 240 hour solutions 
show the circulation continues to move toward the south. Notice again, that the 
circulation does not move too far north. The overall magnitude of the vectors is 
unrealistic at the end of the time period. The large magnitude is caused by the 
arbitrary amplitude factor used in the forcing function. The results of this 
experiment are more consistent with the midlatitude beta-plane theory than the 


results of Experiments 1 through 3. 


«s Z, Z 
о © © 
С) Се) 2 o 
Q 0 со 
СОСКЕ сае эса Чаага 
"иа У 
4 
ЕА 
Пе 


~ 
eee Sax Nip re sese 


NS | 
= Хы 


WES . 
жес RD SEN 
g ESSE 
wee 7. » D i 
КЕ SSS 
LN 


АЛАТ 


p 


11۹ 


За A m 


“2” 


AA 
as oN 
^ 55 


TN 
eevee СЕ 
НДР Иа 


Í 


У т. 


SO www. 7n e @ese ene 
SYS | 


»* 95646 чш MNS 


Fw o... к 
9 99...., ааа ` 


AI 


“2222 | 
е МЕ “а 
Cs fd i ду x 
sy ae zirr SS 
js PIRE SESS 
AAA —— 


аһ 


S 


RSS 
он, "7 у хе 


| 


Wind Vectors at 24 Hours for Phase I, Experiment 4 
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Wind Vectors at 48 Hours for Phase I, Experimen 
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Wind Vectors at 72 Hours for Phase I, Experiment 4 


Figure 6.27. 
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Wind Vectors at 96 Hours for Phase I, Experiment 4 


Figure 6.28. 
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Wind Vectors at 120 Hours for Phase I, Experiment 4 
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Wind Vectors at 144 Hours for Phase I, 


Figure 6.30. 
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Figure 6.32. Wind Vectors at 192 Hours for Phase I, Experiment 4 
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The overall magnitude of the wind vectors of the barotropic mode at 24 
hours (Fig. 6.35) is five times larger than the magnitude of the wind vectors of the 
baroclinic mode at the same time (Fig. 6.36). The vorticity source is much more 


efficient at forcing the barotropic mode. 
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Figure 6.35. Wind Vectors at 24 Hours for Phase J, Experiment 4, 
Mode 1 
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Figure 6.36. Wind Vectors at 24 Hours for Phase I, Experiment 4, 
Mode 2 


5. Experiment 5 

This experiment is exactly the same as Experiment 4, except the sign of the 
frequency (€) in Eq. 2.12 is changed. With a positive frequency (®), U > O in 
Eq. 2.26 and the circulation should be trapped near the source. The wind velocity 
vectors for the lowest level (o = .938) at 24 and 240 hours are displayed in 
Figs. 6.37 and 6.38. The overall magnitude of the vectors is significantly smaller 
than those of Experiment 4. However, the circulation did propagate past the 
equator, but to a much less extent. The propagation 1s most likely due to fast 
moving gravity waves. The modal analysis is consistent with the complete solution, 


and the results of Experiment 4. 
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Figure 6.37. Wind Vectors at 24 Hours for Phase I, Experiment 5 
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Wind Vectors at 240 Hours for Phase I, Experiment 5 
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6. Experiment 6 

In this experiment the forcing is only added to the vorticity equation 
(Eq. 3.1). The equivalent height (H) is 209.09 m, which is the height of the first 
baroclinic mode. The frequency (w) is 26.13 x 10-7 s-! and the period is 27.8 days. 
The wind velocity vectors for the lowest level (o = .938) at 24 and 240 hours are 
displayed in Figs. 6.39 and 6.40. The circulation patterns show a significant tilt 
toward the southwest, but the circulation does not propagate significantly into the 
Southern Hemisphere. The overall magnitude of the barotropic mode vectors (not 
shown) is again much larger than the baroclinic mode vectors (not shown) at the 


24-hour point. 
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Figure 6.39. Wind Vectors at 24 Hours for Phase I, Experiment 6 
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Figure 6.40. Wind Vectors at 240 Hours for Phase I, Experiment 6 
7. Experiment 7 
In this experiment the forcing function is added to the vorticity equation 
(Eq. 3.1). The equivalent height (H) in Eq. 2.12 is 209.09 m, which is exactly 
equal to the first baroclinic mode. The variable o equals x / W, where W is half the 
width of the forcing function in degrees latitude. This value of a is a singular point 
in the complete propagating solution (Eqs. 2.17 and 2.18). The wind velocity 


vectors for the lowest level (6 — .938) are displayed in Figs. 6.41 to 6.50. 
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Wind Vectors at 24 Hours for Phase I, Experiment 7 


W 


120° W 90° 


Figure 6.41. 
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Wind Vectors at 48 Hours for Phase I, Experiment 7 
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Figure 6.42. 
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43. Wind Vectors at 72 Hours for Phase I, Experiment 7 
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Figure 6.44. Wind Vectors at 96 Hours for Phase I, Experiment 7 
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Wind Vectors at 120 Hours for Phase I, Experiment 7 


Figure 6.45 
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Wind Vectors at 144 Hours 


Figure 6.46. 
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Wind Vectors at 168 Hours for Phase I, Experiment 7 


Figure 6.47. 





Wind Vectors at 192 Hours for Phase I, Experiment 7 
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Figure 6.48. 





Wind Vectors at 216 Hours for Phase I, Experiment 7 


Figure 6.49. 
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Wind Vectors at 240 Hours for Phase I, Experiment 7 


Figure 6.50. 


77 


The overall magnitude of the wind vectors 1s approximately the same as in 
Experiment 6, but the circulation pattern has more of a northeast-southwest tilt. 
The frequency (®) of this forcing function is —3.05 x 10-7 s-! and the period is 
238.6 days. The westward movement of the circulation pattern is much slower than 
in Experiment 6, due to the longer period. The modal analysis is consistent with the 
analyses of Experiments 4 through 6. 

8. Phase I Summary 

The results of Phase I indicate that forced Rossby wave forcings applied to 
the thermal equation do not produce significant Rossby-type cold surge responses. 
However, when the same forcings are applied to the vorticity equation, the 
responses are consistent with Rossby-beta-plane theory. Forcing the thermal 
equation consistently produces a stronger response in the baroclinic mode, and 
forcing the vorticity equation produces a stronger barotropic response. When the 
sign of the frequency (w) is changed, the solution does propagate southward. The 
southward propagating solution is not consistent with Rossby wave theory. The 
propagating waves in this case could be gravity waves since the overall magnitude 


of the velocity vectors 1s small. 


C. PHASE II RESULTS 

The forcing functions of Phase II use the impulse function described in 
Chapter V. The maximum amplitude of the forcing is at 24 hours and 6 = .375. In 
each experiment the forcing is applied between 24° N and 36° N, and the model is 
again integrated forward in time for 240 hours. The maximum amplitude of the 
forcing occurs at 24 hours (as shown in Fig. 4.3). The two experiments of Phase II 
are summarized in Table 6.5. Experiments 8 and 9 will compare the effects of the 


impulse forcing in either the thermal or vorticity equations. 
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TABLE 6.5. PHASE II EXPERIMENTS 
Experiment Thermal Forcing Vorticity Forcing 


8 Х 
9 Х 


1. Experiment 8 





In this experiment the impulse source is only added to the thermal equation 
(Еа. 3.4). The wind velocity vector fields for the lowest level (o = .938) are 
displayed in Figs. 6.51 through 6.60. After 24 hours a wavenumber three 
circulation develops near the source latitudes. The overall magnitude of the 


velocity vectors is relatively small and the circulation does not propagate very far 





south. At 48 hours the magnitude increases and the circulation pattern moves 
toward the west, but the circulation is still confined to the source latitudes. By 144 
hours there is a slight northeasterly flow toward the equator. The magnitude of the 


northeasterly flow does not significantly increase after 240 hours. 





Figure 6.51. Wind Vectors at 24 Hours for Phase II, Experiment 8 


19 





Figure 6.53. Wind Vectors at 72 Hours for Phase II, Experiment 8 
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Fi 
igure 6.55. Wind Vectors at 120 Hours for Phase II, Experiment 8 
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Figure 6.57. Wind Vectors at 168 Hours for Phase II, Experiment 8 
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Figure 6.59. Wind Vectors at 216 Hours for Phase II, Experiment 8 
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Figure 6.60. Wind Vectors at 240 Hours for Phase II, Experiment 8 


The overall magnitude of the velocity vectors of the barotropic mode (not 
shown) is on the order of 0.5 m-s-!. The overall magnitude of the first baroclinic 
mode velocity vectors (not shown) is ten times greater. The larger baroclinic 
response 1s consistent with Phase I. The circulation patterns of the modal analyses 
are not significantly different than the total patterns shown in Figs. 6.51 through 
6.60. 

2. Experiment 9 

In this experiment the impulse function is only added to the vorticity 
equation (Eq. 3.1). The wind velocity vector fields for the lowest level (o = .938) 
are displayed in Figs. 6.61 through 6.70. After 24 hours a wavenumber three 
circulation develops and significant circulation patterns propagate toward the 
equator. By 48 hours the magnitude doubles and the northeasterly flow intensifies. 


By 72 hours there is significant flow past the equator. A definite northeast- 
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southwest tilt in the circulation pattern is evident at 96 hours. This tilt is much 
more pronounced than in Experiment 8. By 192 hours there are closed circulation 
cells in the tropics. It is interesting to note that the circulation patterns of this 


experiment also propagate north of the forcing region. 
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Figure 6.61. Wind Vectors at 24 Hours for Phase II, Experiment 9 
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Wind Vectors at 48 Hours for Phase II, Experiment 9 


Figure 6.62 
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Figure 6.63 
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Wind Vectors at 96 Hours for Phase II, Experiment 9 


Figure 6.64. 
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Wind Vectors at 120 Hours for Phase II, Experiment 9 
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Figure 6.65. 
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Figure 6.66. Wind Vectors at 144 Hours for Phase II, Experiment 9 


Len 
e can. MIL ANW 


шу 


ү" 
ЛЫ ^ \\\х+ 
уе h S 
D Jmm 

CUI Nam 

МИН VIAL 
| ИША 
woo 
AN 


We Ж 
حلا‎ ж? fe, 


* e э © чылыы No 


...., да 


44441111. 


„оо AA 


ete 


000444 4 


“ЦГ | 
: iiv 


...7 014. E NN 4 
` 7, | 
соо ЛАК /1 ЖҮ ANANA 


ffi \ 
о 9 9 Зе . Af etn 
Сй Ae 


Ке 
alh; Nn 
ulti m 


4 Î 
ath, Nc 
h QS 
"NM 


š pr 12 
эе уз» J sak \\ 
б ос 


| qM 


=>. 


Ра. 


дал 929 
Nieto? 


oe АА АА, 


“2.444 0005.99 ° 


.% eee Qe 
| ow 


014 


| .. ла 


РФ зо о 


.... А" f 


„Ср рое о оь 


. “өте 
e.. At 


asel}??? 
} 
er mens 


M ss “%%%.7.. 


ео ро... 


4...” 


"9.94 ae 
c" Ме. 


КОРУУ, PRevee., 


5 


60° w 


W 


°w 9 


120 





t 9 


imen 


Wind Vectors at 168 Hours for Phase II, Exper 
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Wind Vectors at 192 Hours for Phase II, Experiment 9 


Figure 6.68. 
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Figure 6.70. Wind Vectors at 240 Hours for Phase II, Experiment 9 
The velocity vector patterns of the first and second modes are similar to 
the complete vector patterns. The overall magnitude of the barotropic mode 
vectors is ten times greater than the baroclinic mode vectors at 24 hours. The 
larger response of the barotropic mode is consistent with the results of Phase I. 
3. Phase II Summary 
The results of Phase II indicate that the thermal impulse forcing is not an 
efficient method of producing a cold surge response. After 240 hours a small but 
significant northeasterly flow develops, and circulation patterns are evident near 
the equator; however, ten days is too slow to be realistic. The impulse forcing 
applied to the vorticity equation does produce a significant, timely response. The 
northeast-south west tilt in the circulation patterns is similar to Bashford (1985) and 


Lim and Chang (1981). The results are consistent with Phase 1, in that the vorticity 


90 


forcing 1s more efficient 1п exciting a barotropic response, and the thermal forcing 


a baroclinic response. 


D. PHASE III RESULTS 

The forcing functions of Phase III use the generic cyclone described in 
Chapter V. In each experiment the baroclinic zone is applied between 30° N and 
60° N. This latitudinal range is representative of a synoptic wave moving through 
the planetary wave over East Asia. The latitudinal structure of both the thermal 
and vorticity forcing functions are shown in Fig. 6.71. The magnitude of each 
forcing equals zero in the Southern Hemisphere and outside the source region. The 
model is integrated forward in time for 240 hours. The maximum amplitude of the 
forcings occur at the 24-hour and 216-hour points (Fig. 6.72). In the vertical, the 
maximum amplitude of both thermal (Fig. 6.73) and vorticity forcings (Fig. 6.74) 
occur at the surface. The three experiments in Phase III will determine which 
forcing is dominant in the generation of cold surges. The experiments are 


summarized in Table 6.6. 


TABLE 6.6. PHASE III EXPERIMENTS 


Experiment Thermal Forcing Vorticity Forcing 
Di X 


X 
X 
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Figure 6.71. Latitudinal Structure of Phase III Forcings 
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Figure 6.72. Temporal Structure of Phase III Forcings 
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Figure 6.73. Normalized Vertical Structure of Thermal Forcing of 
Phase III 
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Figure 6.74. Normalized Vertical Structure of Vorticity Forcing of 
Phase III 
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1. Experiment 10 

In this experiment the forcing 1s applied to both the thermal (Eq. 3.4) and 
vorticity (Eq. 3.1) equations. The wind velocity vector fields for the lowest level 
(6 = .938) are displayed in Figs. 6.75 through 6.84. After 24 hours a wavenumber 
three circulation develops near the source latitudes and there is some propagation 
toward the south. By 48 hours a five meters per second northeasterly wind, typical 
of an East Asian cold surge, is evident to approximately 10° N. After 72 hours the 
northeasterly flow extends past the equator into the Southern Hemisphere. The 
northeast-southwest tilt in the circulation south of the source latitudes, and 
northwest-southeast tilt north of the source is evident at this time. The tilted 
circulation pattern is characteristic of propagating Rossby waves. The propagating 
wave continues in this manner, and at 240 hours there is identifiable circulation at 
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Figure 6.75. Wind Vectors at 24 Hours for Phase III, Experiment 10 
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Figure 6.76. Wind Vectors at 48 Hours for Phase III, Experiment 10 
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Figure 6.77. Wind Vectors at 72 Hours for Phase III, Experiment 10 


95 


60° E 


30° E 
= 2010, 





ге Qao 
eie 
SM diua .. 
M. Стеж 
а ава 


deu 


ine r^ Ж” ј ..і44а» 


WT 
17 


Wind Vectors at 96 Hours for Phase III, Experiment 10 


Figure 6.78. 
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Wind Vectors at 120 Hours for Phase III, 


Figure 6.79. 
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Figure 6.80. Wind Vectors at 144 Hours for Phase III, Experiment 10 
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Figure 6.81. Wind Vectors at 168 Hours for Phase III, Experiment 10 
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Wind Vectors at 192 Hours for Phase III, Experiment 10 


Figure 6.82. 


№ | |, Bu 
\\ Jt MUT 


xu 


ах. 
е 


28 о! ! хум 


Va Wo и 
PEST 


е 
чеди: 


ee Te ие и 2 
KANAA... 
QNNM 


E" ae sce .-тегеттеь 
Mites. ced .....ç..... 


Е 
Е 


NC O 


......,,,/4......... 


| 


چ سو ی 
pice ese = "529-97‏ 


Ps 


E 
| 
| 
| 


4 
via ese esos... 


| 
нелер сеа 


а --5-24262%е0 а о 
4 


ea 
Ее. әӨетеезее е 


irs emet 


b er сс» 


porum 


ISAs свое 


етее онор Дорообо ооо е 


1444%%4»4 1 ee 


} 


еее ооо оета 


Е 


О ЕНЕ осе 


боого љо Ф 2 


тоавачео 


хосте 


Фтеғе.ье 


@etesean 


вооот оо 


Тато ае 


e=“. .... 


про зе 


гето ава 


eeeet se 


ue oe FO 


....... 


*e228242209 


**^.* ve. 


@@c =e... 


v^. ^on. 


зоо 





{ 10 


пеп 


Ехрег 


Wind Vectors at 216 Hours for Phase III, 


Figure 6.83. 


98 


aes 


Е. 


Ne 
N es 
= 
appa 


|| | 
е Како 


eaP Angee 
4... Vie E +00440 
+! 


1 
veter — — си | > ль. + * 


i 


Це 
EU 

“A ANS NS 

...4...- See пи А 
= ал 34311 one а > TENUES f) 


* NA er^ ~ 
' 
' 


i 
` 
2 
^ 
4 
, 
е 
РА 
r 
° 
+ 


Ах 


.....4.%%, 


| 
ме реа 


i 
ace eda АҚ — 


| 
90° W 60° W 





Figure 6.84. Wind Vectors at 240 Hours for Phase III, Experiment 10 

The barotropic (Fig. 6.85) and baroclinic (Fig. 6.86) mode analyses at 24 
hours show similar circulation patterns as the complete solution at 24 hours. The 
overall magnitude of the barotropic mode vector field is five times larger than the 
baroclinic vector field. The larger barotropic response suggests the vorticity 
forcing is dominant in the generic cyclone. The higher modes have similar 


circulation patterns, but the overall magnitude of the vector field is much smaller. 
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Wind Vectors at 24 Hours for Phase III, Experiment 10, 


Figure 6.85. 


Mode 1 





Wind Vectors at 24 Hours for Phase III, Experiment 10, 


Figure 6.86. 


Mode 2 
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2. Experiment 11 


In this experiment the forcing is only applied to the thermal equation 
(Eq. 3.4). After 24 hours a wavenumber three circulation develops in the lowest 
level (o = .938), but the circulation is rather chaotic (Fig. 6.87). The circulation 
remains trapped for the whole time period. After 240 hours (Fig. 6.88) the 
sinusoidal circulation pattern is still within the source latitudes. It is apparent that 
the thermal forcing does not play a significant role in producing cold surges in this 
model. The overall magnitude of the baroclinic mode vector field at 24 hours (not 
Shown) is four times larger than the magnitude of the barotropic field. The larger 


baroclinic response is similar to the results of Phase I and II. 





Figure 6.87. Wind Vectors at 24 Hours for Phase III, Experiment 11 
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Figure 6.88. Wind Vectors at 240 Hours for Phase III, Experiment 11 


3. Experiment 12 
In this experiment the forcing is only applied to the vorticity equation 
(Eq. 3.1). Since the vorticity forcing is dominant, and the thermal forcing does 
not propagate away from the source, the circulation patterns and magnitudes are 
not significantly different than the combined thermal and vorticity forcing of 
Experiment 11. 
4. Phase III Summary 
The results of Phase III indicate the vorticity forcing is the dominant 
component in the nonlinear interaction of a planetary wave and a synoptic-scale 
generic cyclone. The thermal-only forcing does not propagate away from the 
source and is less intense than the vorticity-only forcing. The vorticity forcing 
(and thus the combined thermal and vorticity forcing) does produce a northeasterly 


wind within 48 hours, which is typical of East Asian cold surges. The northeasterly 
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flow continues into the Southern Hemisphere, which suggests a tropical response to 


a midlatitude forcing (Lim and Chang, 1987). 


103 


VII. CONCLUSION 


The objective of this study was to determine whether the nonlinear interaction 
between a planetary-scale wave (wavenumber four) and a rapidly growing 
synoptic-scale wave (wavenumber seven) could force a northeasterly cold surge. 
The nonlinear interaction of these waves produced a wavenumber three forcing. 
The basic model was a linearized global spectral model with eight levels. Fourteen 
experiments with three different types of source functions were applied to the 
thermal and/or vorticity equations. The solutions of each experiment were 
projected onto the vertical normal modes of the primitive equation model. 

The results of each phase of the experiments indicate a midlatitude thermal 
forcing does not produce a significant cold surge response. The impulse forcing, 
similar to Lim and Chang (1981), applied to the thermal equation produced a weak 
cold surge response after 240 hours; however, ten days is longer than the typical 
two-day surge generation period observed in the East Asian winter monsoon. The 
difference between the equatorial beta-plane model used by Lim and Chang (1981), 
and the spherical primitive equation model used in this study will require further 
explanation. However, the thermal forcing consistently produced a stronger 
response in the baroclinic mode. 

The midlatitude vorticity forcing, on the other hand, produced a stronger 
response in the barotropic mode. The wavenumber three vorticity forcing 
produced a significant cold surge response within 48 hours. The cold surge 
produced by the vorticity forcing is consistent with the observed cold surge in the 


East Asian winter monsoon. 
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The lack of cold surge response to a midlatitude thermal forcing, and the 
significant response to a midlatitude vorticity forcing may be explained by 
geostrophic adjustment theory. In midlatitude geostrophic adjustment the wind 
field tends to adjust to the mass field. Thus, when the forcing is prescribed by 
heating, the local adjustment (to rotational motion) is more complete. When the 
forcing is prescribed by vorticity, the adjustment will be less efficient and a 
considerable amount of energy will propagate away in divergent motion. The 
equatorward propagation studied by Lim and Chang's (1981) theory is through 
equatorial Rossby waves that have a significant divergent component. Thus, the 
divergent motion may offer a possible explanation for the different responses to 
the two types of forcings. 

The initial condition of this study did not include any basic current (mean 
flow). In reality, the cold surge is superimposed on an existing northeasterly 
monsoonal flow at the surface. It would be interesting to see if the presence of an 
initial mean flow would alter the results of the forcings. Also, the nonlinear effects 
of topographic channeling play an important role in the cold surge response. Using 


some type of topographic barrier may lead to interesting conclusions. 
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